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Abstract. High resolution OVRO aperture synthesis 
maps of the ^^CO 1-0 emission in the "Medusa" galaxy 
merger (NGC 4194) reveal the molecular emission being 
surprisingly extended. The ^^CO emission is distributed 
on a total scale of 25" (4.7 kpc) — despite the appar- 
ent advanced stage of the merger. The complex, strik- 
ing ^^CO morphology occupies mainly the center and the 
north-eastern part of the main optical body. The extended 
^^CO flux is tracing two prominent dust lanes: one which 
is crossing the central region at right angle (with respect 
to the optical major axis) and a second which curves to 
the north-east and then into the beginning of the northern 
tidal tail. 

The bulk of the ^^CO emission (67%) can be found in a 
complex starburst region encompassing the central 2 kpc. 
The molecular gas is distributed in five major emission 
regions of typical size 300 pc. About 15% of the total ^^CO 
flux is found in a bright region I'.'S south of the radio 
continuum nucleus. We suggest that this region together 
with the kpc sized central starburst is being fueled by gas 
flows along the central dust lane. We discuss the merger 
history of NGC 4194 and suggest that it may be the result 
of a early- type/spiral merger with a shell emerging to the 
south of the center. 

The total molecular mass in the system is estimated 
to be at most 2 x 10^ M©, depending on which ^^CO - 
H2 conversion factor is applicable. The high ^-^CO/^^CO 
1-0 intensity ratio, ~ 20, indicates highly excited physical 
conditions in the interstellar medium showing that the 
starburst has a big impact on its surrounding ISM. At the 
current rate of star formation, the gas will be consumed 
within 40 million years. 
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1. Introduction 

The Medusa merger, NGC 4194, (Table 1; Fig. 1) be- 
longs to a class of lower luminosity (Lfir — 8.5 x 
10^" L0 at D=39 Mpc) mergers, an order of magni- 
tude fainter than well known Ultra Luminous Galaxies 
(ULIRGs (Lfir- 10^^ Lq)) such as Arp 220, but with an 
extended region of intense star formation (Prestwhich et 
al. 1994; Armus et al. 1990). As the name indicates, the 
morphology of the Medusa merger is as spectacular as 
that of the brighter objects. The "hair" or "tentacles" are 
a patchy tidal tail stretching out 60" (11.3 kpc) north of 
the main body, which appears twisted. There also is a 
sharp, curved feature (and a very faint tail) to the south. 

NGC 4194 has been mapped in ^^CO 1-0 and 2-1 with 
a single dish telescope by Casoli et al. (1992) and de- 
tected in the more highly excited ^^CO 3-2 line by Dev- 
ereux et al. (1994). It is one of the handful known mergers 
(e.g. Arp 220, IC 694, NGC 6240) with unusually high 
12QO/13CO 1-0 intensity ratios, 7^ >20 (e.g. Aalto et al. 
1991b; Casoli et al. 1992). In contrast, the majority of 
large galaxies exhibit TZ in the range 5 — 15 (e.g. Aalto 
et al. 1995). One feature these HR (High Ratio) mergers 
have in common is a high FIR /(60/xm)//(100/ini) flux 
ratio iiO.S - which likely is an indication of a high dust 
temperature (e.g. Aalto et al. 1991a, 1991b). High resolu- 
tion ^^CO maps of several luminous (Lir>10^^ Lq) HR 
mergers reveal that nearly all the molecular gas is con- 
centrated within a few hundred parsecs of the center (e.g. 
Scoville et al. 1989; Bryant & Scoville 1996; Downes & 
Solomon 1998). 

The resulting large central surface densities require 
high pressures in hydrostatic equilibrium and the molec- 
ular ISM becomes highly turbulent. Large turbulent line 
widths in molecular clouds will result in reduced optical 
depths in the ^^CO 1-0 line and the i^CO/i^CO 1-0 line 
ratio will be elevated compared to the optically thick ratio 
(for a discussion see Aalto et al. 1995). So for the lumi- 
nous HR mergers this could be one important mechanism 
behind the large observed ^^CO/^'^CO line ratios (Aalto 
et al. 1995), aside from the additional effect from large ki- 
netic temperatures. Will the lower luminosity HR galaxy 
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NGC 4194 also follow this pattern of a compact nuclear 
starburst, or will an extended ^^CO distribution require a 
different dominant mechanism behind the elevated TZ? 

Most of the high resolution imaging effort of -'^^CO in 
mergers has been directed toward the ULIRGs and not 
much is known, statistically, of the more modest luminos- 
ity merging galaxies. Apart from addressing the question 
of the compactness of the ^^CO emission in an HR merger, 
a ^^CO image of NGC 4194 will improve the statistics of 
the "'^^CO distribution in moderate luminosity mergers. 

The elevated values of TZ in hot dust mergers have also 
been interpreted as a global underabundance in ^^CO with 
respect to ^^CO and not an effect of reduced optical depths 
in ^^CO (e.g. Casoli et al. 1992). However, the moderate 
optical depth (r w 1) in ^^CO 1-0 towards centers of 
galaxies in general, and ULIRGs in particular, is now a 
reasonably well established fact (e.g. Aalto et al. 1997; 
Dahmen et al. 1998; Downes & Solomon 1998). 

We have imaged NGC 4194 in the ^^CQ 1 line with 
OVRO to study the distribution and kinematics of the 
molecular gas in an HR merger galaxy of moderate lumi- 
nosity. We have also obtained single dish spectra of the 
1-0 transitions of ^^CQ, i3C0 and HCN. The observa- 
tions are presented in Sect. 2; the results in Sect. 3; the 
consequences of these findings are discussed in terms of 
merger history, central gas surface densities and final fate 
of the molecular gas in Sect. 4. 

2. Observations 

We have obtained high resolution maps of ^^CO 1-0 using 
the Caltech six-element Owens Valley Radio Observatory 
(OVRO) millimeter array. Two tracks were taken in the 
low resolution configuration in September and October 
1996 and two in the high resolution configuration in De- 
cember 1996. The naturally weighted synthesized beam 
size is 2 '.'5 x 2'.'0 (472 x 378 pc for D=39 Mpc) and the 
beam position angle (BPA) is -66°. The sensitivity of this 
map is 11 mJy beam~^ channel"^, corresponding to 0.19 
K channel"-'^. For a beamsize of 2 "3 a brightness temper- 
ature (Tb) of 1 K corresponds to 57 mJy at A=2.6 mm. 

For the (naturally weighted) high resolution data only, 
the beamsize is 2'.'0 x 1'.'7 and BPA is -71°. For a resolu- 
tion of 1 '.'85 a brightness temperature of 1 K corresponds 
to 36 mJy. 

System temperatures were 500 600 K. The quasar 
1150-1-497 was used for phase calibration and Neptune and 
Uranus for absolute flux calibration. The digital correla- 
tor, centered at Vlsr = 2560 kms""'^, provided a total 
velocity coverage of 1170 kms^^. Data were binned to 4 
MHz resolution, corresponding to 10 kms~^, and to con- 
struct the map, the resolution was reduced to 20 kms~^. 

Single dish spectra of the 1-0 transition of ^^CO, ^'^CO 
and HCN were obtained toward the center of NGC 4194 
with the OSO 20m telescope. The observations were car- 
ried out in October 1999, and the system temperatures 
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Fig. 1. Contour plot of an R-band CCD image of 

NGC 4194 (Mazzarella & Boroson 1993). The levels arc 
logarithmic to clearly show the faint structure, such as the 
northern tidal tail. 



Table 2. 


Sin^ 


^le dish 


results 


Line i 


H") 


/ Tmb 


(K kms-i) 


12 CO 


33 




17 ±2 




34 




0.9 ±0.2'' 


HCN 


43 




<0.4'' 



a) : The i^CQ/i^cO 1-0 intensity ratio is 19 ± 4. 

b) : The i^CO/HCN 1-0 intensity ratio is > 25, corrected for 
beamsize. 

were typically 400 K for ^^CO and ^^CO and 200-300 
K for HCN. The pointing accuracy was checked on SiO 
masers and was found to be about 3". Beam sizes and 
efficiencies are listed in Table 2, together with the results. 

3. Results 

3.1. Morpholology and kinematics of the molecular 
distribution 

Integrated intensity: The naturally weighted "'^^CO emis- 
sion (Fig. 2a) is distributed on a scale of 30" (5.7 kpc). 
The bulk of the emission (67 %) is concentrated in an 10" 
(1.9 kpc) sized oval structure centered on the optically 
bright main body. Fainter, patchy emission curves out to 
the north from the brighter emission. 

An overlay of the ^^CO distribution on an R-band im- 
age (Fig. 3) reveals that the general ^^CO morphology is 
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Table 1. Properties of NGC 4194^) 
R.A. (1950.0) Dec (1950.0) Type rub Size (blue) PA Lfir 



/(60/jm) 
/(lOOMm) 



M(HI) 



M, 



dyn 



12'^ 11"" 41.22= 54° 48' 16'"'' Pcc'=' IS.O'^' 2'.3 x l'.6 146' 



d) 



0.97 



2.2 X 10'-' Mg 7 X 10^ 



a) : Other aliases: UGC 07241, Arp 160, Mrk 201, VV 261, I Zw 033 

b) : Radio continuum peak position (Condon et al. (1990) 

c) : from the UGC catalogue (Nilson 1973) 

d) : as measured the bright main body of the optical R-band image (from Mazzarella & Boroson) 

e) : (lO/im - lOO^m), Dcutsch & Willner (1987) (for D=39 Mpc) 

f) : Thuan and Martin (1981) (for D=39 Mpc) 

g) : Adopting the inclination i « 40° and using Fig. 4, wc get Kot=115 kms~^. The dynamical mass is 
Afdyn ~ 2.3 X 10*(-^^)'^( j^) M0 where Viot is in kms^^ and R in pc. For the inner 8" (1.5 kpc) radius the dynami- 
cal mciss is Mdyn = 7 X 10^ M©. In this unrelated system, one should of course use dynamical masses, derived from the 
assumption of regular rotation, with care. 
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Fig. 2. a) Integrated intensity of the naturally weighted ^^CO map. Contour range is (0.54, 1.6, 2.7, 3.8, 31) Jy 
beam^^ kms^^. The grayscale ranges from 0.7 (light) to 27 (dark) Jy beam^^ kms^^. The radio continuum position 
is marked with a cross. The total flux in map is 170 Jy kms^^. The dashed line indicates the location of the position 
velocity cut in Fig. 6 b): The velocity field. Contours range from 2400 to 2640 kms~^, with spacing 24 kms~^. The 
grayscale ranges from 2400 (light) to 2650 (dark) kms~^ 



quite different from that of the main optical body of the 
galaxy. The molecular gas extends from the center to the 
north-eastern part of the object, and the major axis of the 
central CO intensity distribution is almost perpendicu- 
lar to that of the optical body. The position angle (PA) of 
the ^^CO oval intensity distribution of the bright emis- 
sion is « 49° and the FWHM source size is 4 '.'2 x 2 '.'5 The 
position angle of the optical emission is 146°. 

The curved northern ^^CO emission reaches into the 

base of the tidal tail, and some faint emission is also de- 
tected near a bright patch (indicated with an arrow in 
Fig. 3) in the tail itself. In Fig. 3, the ^^CO emission ap- 
pears to be inside of an optical protrusion, or arm, in the 
east (marked with an arrow in Fig. 3). Armus et al. (1990) 



describe this feature as a "spur" or a "jet" . The blue band 
image in Arp's Catalog of Peculiar Galaxies (1966) gives 
instead a strong indication of a foreground eastern dust 
lane blocking the background emission and creating the 
impression of a protruding feature separate from the main 
body. The molecular emission traces the space between 
the "spur" and the main body which is consistent with 
the interpretation of a dust lane. 

The total ^^CO flux in NGC 4194 detected by OVRO 
is 170 Jy kms""'^ which, for a standard "'^^CO to H2 mass 
conversion factor, translates to 2.3 x 10^ Mq (0.92 x 
IQ^SAVD"^ M© for X = N(H2)/I(i2co) = 2.3 x lO^o 
cm~^ (K kms"^)""'^) of molecular gas. Within calibra- 
tional error bars, we recover all of the single dish flux. 
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Fig. 3. Overlay of the ^^CO contours (white, apart from 

clouds in the tail marked in black) on a greyscale, overex- 
posed optical R-band image (Mazzarella & Boroson 1993). 
The left arrow indicates the eastern expansion and the 
right one points to one of the bright condensations in the 
northern tidal tail. 

The peak flux is 27 Jy beam~^ kms~^ (16% of the total 
flux) at a: 12'' ll'" 41.22" ; 5: 54° 48' 14'.'4. which trans- 
lates to a molecular mass of 3.8 x 10^ M0 within a radius 
of 240 pc. However, we suspect that the conversion fac- 
tor overestimates the molecular mass in this system — in 
particular in the center (see Sect. 4.3.1.). 

The velocity field: In Fig. 2b we present the ^^CO veloc- 
ity field which shows the complex dynamics of NGC 4194. 
Close to the center there seems to be organized rotation 
of the gas (even though the velocity contours are far from 
a regular spider diagram) with a position angle of 160° — 
close to the optical major axis. Thus, the major axis of the 
rotation of the gas follows the PA of the optical light distri- 
bution, rather than the distribution of the ^^CO emission. 
A position- velocity (pV) cut at PA 160° (Fig. 4) indicates 
a solid-body rotation curve which starts to flatten (or fall) 
at a radius of 2" (380 pc). The Figure also shows that the 
^^CO peak is offset from the kinematic center by 1". The 
projected rotational velocity is 90 kms"'^. Thus, for an 
inclination of 40° (from the optical isophotes) the depro- 
jected rotational velocity is 140 kms""' in the inner 16" 
(3 kpc). The systemic velocity at the kinematic center is 
2510 kms"\ 

The extended northern gas in general does not move 
according to the pattern established in the inner region. 
The eastern part of the emission is at higher velocities 
than the western emission and velocity discontinuities can 
be traced along the curved ridge at the transition between 
the velocity systems. The highest velocity in the system, 
2700 kms""', can be found on the eastern dust lane at 5: 
54° 48' 22". The velocity of the gas in the eastern dust 
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Fig. 4. pV cut at PA 160° across the central ^^CO body. 

The location of the cut is indicated in Fig. 2a. The cross 
marks the estimated kinematic center of a solid body ro- 
tating structure. 

lane then gradually drops along the lane, until it reaches 
the tidal tail where the velocity of the gas clumps is 2580 
kms"'^. 

The Dispersion Map: Fig. 5 displays the dispersion 
map with a few selected spectra. The one dimensional dis- 
persion is quite moderate, aidS 10 kms^^ in large parts 
of the galaxy. Regions of high dispersion are found where 
the velocity field shows discontinuities or irregularities, 
and the sample spectra reveal that here the enhanced ve- 
locity dispersion is caused by multiple spectral features 
rather than genuinely broad lines. The velocity shifts be- 
tween the narrow peaks is 150-200 kms^^, i.e. the typical 
rotational velocity of a spiral galaxy. There is also a con- 
tinuous curved region of elevated dispersion (i.e. multiple 
spectra) in the center of the map. Even the width of the 
spectrum at the intensity peak is quite low (lower panel, 
right hand side) with (Tid = 50 kms""', which is narrow 
for the center of a merging system. 

Channel Maps: Dividing the emission into three ve- 
locity bins, (2380 2485 kms"^), (2506 - 2570 kms"^), 
(2590 - 2720 kms^^), (Figs. 6a,b,c) shows the emission be- 
ing distributed in long, arm-like features of quite different 
appearence in the three bins. The low and high velocity 
bins were chosen to include each one of the two features 
seen in the double peaked spectra (e.g. middle right panel 
in Fig. 5). The middle bin includes emission at intermedi- 
ate velocities near the systemic value. 

The low velocity map (Fig. 6a) shows a patchy distri- 
bution extending over 20" to the north-south. Closer to 
the center, brighter emission curves 5" in the east-west 
direction. 

At mid- velocities (Fig. 6b), the distribution is dom- 
inated by a striking, elongated curved shape, approxi- 
mately 15" (2.8 kpc) in length. This structure is the main 
cause of the discrepancy between the ^^CO emission dis- 
tribution and the R-band light distribution. It traces a 
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Fig. 5. Dispersion map witii sample spectra. The contours are (12, 36, 60)kms ^ and the grayscale ranges from 6 



(light) to 60 (dark) kms" 
central region. 



The sample spectra are selected to show regions of multiple spectral features and the 



central prominent dust lane. To show this, we have over- 
layed the mid- velocity ^^CO emission on an archival HST 
image in Fig. 7. Another, fainter, 20" (3.8 kpc) long curved 
structure is found in the north. This feature may be a 
continuation of the central dust lane and the gap between 
them a depression in the '^^CO brightness. 

The high velocity emission bin (Fig. 6c) shows a bright 
peak in the center with a source size of 2 '.'5 x 1 '.'8. A promi- 
nent, narrow curved feature stretches out to the north- 
east, (20" to the north) and follows most of the eastern 
obscuring dust lane seen in Arp's B-band image. We sug- 
gest the northern, curved feature in the mid-velocity bin 
is the continuation of this curve. As mentioned above, the 
central dust lane may connect to this feature as well. In 
Fig. 6d the features in the three velocity bins are marked 
with dashed lines. 




3.1.1. The High Resolution Map 
center 



a close-up on the 



For a closer look at the structure of emission in the inner 
10", we selected the high resolution data only, and thus 
increased the resolution to 1'.'7 (but lost some sensitivity 
and 7iu-covcragc). The integrated intensity map (Fig. 8a) 
of the inner region consists of 5 peaks, a, b, c, d and e (Ta- 
ble 3) and extended emission surrounding them. The total 
flux of the inner 10" x 10" is 64 Jy km/s. The brightest 
peak, a (17 Jy kms~^ beam~^), is at the same position 
as in the lower resolution map, i.e. close to the radio con- 
tinuum peak, but shifted to the south by 1"5 (280 pc). 
Feature a has a size of I'.'S x 1'.'3 and a position angle 



Fig. 7. Contours of the mid- velocity (2506 2570 km s^ ) 
"^^CO emission overlaid on an archival HST WFPC2 im- 
age which shows the ^^CO emission tracing a prominent 
central dust lane. 



of 45°. The fitted FWHM sizes of features a-e are typi- 
cally 200-400 pc, which is comparable to the sizes of Giant 
Molecular Associations (GMAs) found in other galaxies. 
Features b and c can be identified with knots in the Ha 
distribution (Armus et al. (1990)). The peak intensity in 
the high resolution map is 200 mJy, which corresponds to 
5.4 K. 
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Fig. 6. The left panel a) shows the integrated intensity (naturally weighted "'^^CO map) in the lower part of the 
velocity range (2380 2485 kms^^), panel b) shows the central part of the velocity range (2506 ~ 2570 kms~^), and 
panel c) shows the highest velocities (2590 - 2720 kms~^). The dashed lines indicate the position of the central dust 
lane from the mid velocity panel. Contours are (0.7, 1.4, 2.1. 2.8 ... 14) Jy beam^^ knis^^. The cross marks the radio 
continuum position. The right panel bf d) is a sketch where the features from the three velocity bins are indicated 
with dashed lines. 




Fig. 8. a) Map of total integrated intensity for the high resolution data only showing the inner region of NGC 4194. 
Contours arc (0.5, 1.5, 3.5, ...) Jy beam^^ kms^^. The peak flux is 17 Jy beam^^ kms^^. The intensity peaks are 
designated a to e with falling peak intensity. The grayscale ranges from 1 (light) to 16 (dark) Jy bcam^^ kms^^. b) 
The high resolution velocity field. Contours and grayscale arc the same as in Fig. 2b. The position of the ^^CO peak 
a is indicated with a black triangle, c) High resolution dispersion map. The position of the ^^CO peak a is indicated 
with a black triangle. 



Fig. 9 is an overlay of the integrated high resolution 
^^CO emission on an archival HST image (enhanced to 
show the dust lane, same as Fig. 5). The peaks a, b, d and 
e all lie on the central dust lane. 

The velocity field (Fig. 8b) is slightly more twisted 
than in the lower resolution map and covers the same ve- 
locity range. Features d and c are causing twists in the 
velocity field because they are at too low velocities to fit a 
spider diagram. The maximum dispersion (100 kms~^) in 
Fig. 8c does not occur in a, but in an elongated (east-west) 
region located north of the ^^CO peak, close to feature b. 
This is caused by the crossing of the dust lane, i.e. by a 
superposition of multiple spectral features. 

Fig. 10 shows an overlay of the ^^CO emission on the 
1.49 GHz radio continuum structure (Condon et al. 1990). 
There is a roiigh agreement between the CO and the 
radio-continuum distribution which seems to break down 



at smaller scales (i.e. less than 1 kpc). The ^^CO peak 
(feature a), for example, is located 1" south of the radio 
continuum peak. There is an arm-like structure curving 
from the center to the west which has a ^^CO counterpart 
in feature c. The one ^^CO feature entirely lacking in the 
continuum distribution is e in the north. Comparisons be- 
tween radio continuum and CO distributions in starburst 
galaxies often show a correlation between the two distribu- 
tions (e.g. Yun et al. 1994, Bajaja et al. 1995). This effect 
is often attributed to cosmic ray heating of the molecular 
gas (e.g. Adler et al. 1991; Allen 1992). We will discuss 
this correlation further in Sect. 4.2. 

3.2. Single dish results 

Results arc presented in Table 2 and ^^CO, ^^CO spectra 
in Fig. 12 (HCN was not detected, and its spectrum is 
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Table 3. Central Emission Peaks 



Peak 


R.A. (1950.0) 


Dec (1950.0) 


Peak Intensity 


Size (FWHM) 


Vc (kms-') 


Notes 




(1950.0) 


(1950.0) 


(Jy kms"'^ beam"'^) 




(kms"') 




a 




54° 48' 14 '.'4 


17 


1'.'5 


2592 


'2 CO peak 


b 


12'' 11"" 41.28' 


54° 48' 15 '.'9 


13 


2 '.'5 X O'.'S 


2531 


on dust lane''^ ; Ha knot 


c 


12'' ll"" 40.82' 


54° 48' 15 '.'4 


8.6 


1'.'7 X 0'.'5 


2496 


Ha knot 


d 


12'' 11"" 40.75'' 


54° 48' 13 '.'9 


6.9 


I'.'S X 0'.'7 


2539 


on dust lane 


e 


12'' 11"" 41.57' 


54° 48' 18 '.'4 


4.6 


2 '.'4 


2563 


on dust lane 



a): Double spectrum: peaks at 2464 km s ^ and 2549 km s ' (main) 
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Fig. 9. Overlay of the integrated high resohition ^^CO 
emission in contours on an archival HST image (same 
as Fig. 5). The central dust lane is indicated as a white, 
dashed hne, and most of the ^^CO peaks (a, b, d, e) lie on 
this lane. The '^^CO peaks are indicated with their corre- 
sponding letters. 

not shown). The ^^CO flux detected by OVRO is, within 
calibrational error bars, the same as measured by the OSO 
20m telescope, which is approximately 60% more than the 
flux Dcvereux et al. (1994) find with the FCRAO 14m 
telescope. 

We note that the integrated ^^CO/^^CO intensity ratio 

we find is lower, « 20, than previously reported by Casoli 
et al. (1992), « 55, although their error bar of ±22 in- 
cludes our measurement. Our result means that NGC 4194 
is not an extreme HR merger and '^^CO 1-0 is quite de- 
tectable, even if the ratio is clearly above what is usually 
found for starburst galaxies. There is a hint of a line shape 
difference in the '^^CO and '^^CO lines with an extra nar- 
row component at lower velocities in the ^'^CO spectrum. 



z 
o 

< 
Z 
J 

o 

td 
O 




08 ^ ^ ^ \ \ 1— 

121142.0 41.8 41.6 41.4 41.2 41.0 40.8 40.6 40.4 
RIGHT ASCENSION (B1950) 



Fig. 10. Overlay of the radio 1.49 GHz radio continuum 
(Condon et al. 1990) in contours on the high resolution 
^^CO emission in grayscale. The ^^CO features are la- 
beled with their corresponding letter. The contours are 
logarthimically spaced, as in Condon et al.. 

The signal to noise ratio in the ^^CO spectrum is too low 
to warrant further speculation on this matter here, but a 
high resolution '^^CO image may reveal interesting excita- 
tion changes across the galaxy, as has already been shown 
for the merger starburst system Arp 299 (Aalto et al. 1997; 
Casoli et al. 1999). The HCN emission is at least as faint 
as the ^^CO emission with the ^^CO/HCN intensity ratio 
5; 25. In comparison to other FIR huninous galaxies (e.g. 
Solomon et al. 1992) the HCN is relatively faint. 

4. Discussion 

4.1. Merger timescale and ^^CO morphology 

4.1.1. The Medusa as a major disk-disk merger 

Joseph & Wright (1985) suggested that NGC 4194 is a 
collision between two disk galaxies and put it fairly late 

in a merger sequence: after Arp 220 and before NGC 3310. 
They used the apparent diffuseness of the tidal tail and 
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Fig. 11. The OSO 20m telescope single dish spectrum of 
^^CO (upper panel) and ^"^CO (lower panel). The intensity 
scale is in Tmb- 



the coalesced optical main body as principal arguments 
for this placement in the evolutionary scheme. 

However, despite the suggested advanced stage of the 
NGC 4194 merger, the molecular gas structure is com- 
plex and distributed on a scale of several kpc. In the late 
stages of a merger, simulations suggest that the bulk of the 
molecular gas would have sunk to the inner 200 pc of the 
merger, regardless of the previous morphology/orientation 
of the precursor galaxies (e.g. Barnes & Hernquist 1996). 
The corresponding size scale for NGC 4194 is an order of 
magnitude greater than that. The funneling of molecular 
gas toward the center of NGC 4194 has been less effective 
than for ULIRGs, such as Arp 220 or NGC 6240, which is 
interesting since both these mergers appear younger than 
NGC 4194. Apart from the ULIRGs being more luminous 
by factors of at least 10, they also have broader line widths, 
and deeper potentials. If NGC 4194 is the result of the 
merger of two disk galaxies of comparable mass, i.e. a 
'major merger', cither the merger is less advanced than 
previously thought, or the intrinsic and/or orbital proper- 
ties of the NGC 4194 progenitors are, after all, important 
factors. 

One can use the presence and separation of a dou- 
ble nucleus as an indication of merger age. In NGC 4194 

there is no obvious sign of another nucleus, but there is 
a fairly bright extension to the south of the radio con- 



tinuum emission, which is the location of the ^^CO peak, 
a. If a is another nucleus in the merging system, then the 
projected separation between them is I'.'S, or 280 pc. How- 
ever, the non-compact appearance (and lower lumniosity) 
of a indicates instead that it is a bright star forming re- 
gion rather than a second nucleus, in which case the upper 
limit of the nuclear separation becomes even more strict. 
It is perhaps also possible that merger orientation hides 
the second nucleus behind the first one, but this should 
be a rather unlikely scenario. 

We find it difficult to reconcile the extended, unre- 
laxed appearance of the molecular gas distribution with 
a major merger scenario where the lack of a double nu- 
cleus indicates a very advanced merger. Together with the 
morphological appearance of the Medusa, this lead us to 
suggest another merger scenario: 

4.1.2. The Medusa as a yoimg shell-galaxy 

The southern sharp, curved feature seen in the optical 
image may well be a shell in its early stages of formation, 
thus raising questions on whether an elliptical/early- type 
galaxy was one of the progenitors. The total molecular 
mass of 2 X 10^ Mq is indeed less than what one would 
expect from the collision of two gas-rich disks. 

Quinn (1984) found that shells can form through 
"phase-wrapping" of dynamically cold disk material in a 
rigid, spherical potential. The similarity between the over- 
all Medusa morphology and simulations of a radial, pla- 
nar collision between a disk galaxy and a more massive 
spherical, is striking (e.g. Quinn 1984, Fig. 1 timestep 6; 
Kojima & Noguchi (KN) 1997, Fig. 4a timestep 4.7, 6.2). 
The first shell is forming on the opposite side of the sin- 
gle, broad tidal tail creating an almost "cometary" struc- 
ture. Although, we note that the implied molecular mass 
of NGC 4194 suggests that the disk galaxy involved was 
more massive than a dwarf — which is often assumed in 
the simulations. 

KN also simulated the behaviour of the gaseous com- 
ponent during the encounter and found that the gas be- 
comes scattered much like the disk stars. This is in sharp 
contrast to the models of Weil & Hernquist (WH)(1993) 
who instead conclude that the gas from the disk galaxy 
quickly sinks to the center of the merger. This also leads to 
different views on the occurence of starbursts in these en- 
counters. In the KN approach, the scattering of gas clouds 
prevents star formation, while WH suggest that starbursts 
are likely to occur in the central gas concentration. In the 
WH model the shells are expected to be gas-poor because 
of the separation between gas and stars. The difference be- 
tween the KN and WH models is largely caused by their 
ISM represenation: KN using sticky particles and WH a 
hydrodynamical model. 

Now, if we adopt the notion that NGC 4194 is indeed a 
radial collision between a spiral and a more massive, gas- 
poor spherical (i.e. an unequal mass, or 'minor' merger). 
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wc can make the following observations: No molecular gas 
is detected in the southern shell in accordance with the 
WH model, but perhaps more sensitive observations can 
find gas here. Both HI and CO has been detected in the 
shells of Centaurus A (Schiminovich et al. 1994; Charman- 
daris, Combes & van der Hulst 2000). The molecular gas 
is not scattered all over the system, nor can all of it be 
found in the very inner region. The distribution is some- 
where in between the two models, and the gas in the inner 
region seems to be going through a starburst on a fairly 
large scale (see Sect. 4. 2). The multiple spectral compo- 
nents may reflect gas on diflferent orbits and in different 
stages of its dance around the center of the elliptical. The 
velocity shift between the two components is roughly the 
same as the typical rotational velocity for a spiral galaxy. 
Thus the existence of two spectral components does not 
necessarily imply that the molecular gas is coming from 
two different progenitor galaxies. For both the major and 
minor merger models, however, the presence of two spec- 
tral components suggests that the gas has some time to 
go before it reaches relaxation. 

However, the elliptical-spiral model leaves us with (at 
least) two quesions: 1) where is the massive elliptical? 
2) what is the significance of the (very) faint tail to the 
south (which can only be seen in deep images)? The first 
question can better be answered once there is a near in- 
frared K-band image that can see through the dust lanes 
and dust enshrouded starburst currently hiding the center. 
The southern tail is very diffuse, and may be a remnant of 
an earlier interaction. Its significance can be better judged 
by sensitive HI observations. 

We have made a search through galaxy morphology 
catalogues looking for Medusa-like mergers. We found 
three more objects with similar morphology (i.e. one 
prominent tidal tail and on the opposite side of the galaxy 
a sharp curved feature): ESQ 240-GOl, NGC 7135 and 
NGC 5018. So, the Medusa morphology is not unique, 
but fairly rare, which could be an effect of the rather short 
timescales involved in the radial encounters. The Medusa 
still stands out among these objects because of its bright 
FIR emission and hot dust. Perhaps the Medusa is in a 
short-lived stage of intense star formation that will con- 
sume a large fraction of the molecular gas, leaving the 
galaxy to look more like NGC 7135 once the burst is over. 

4-2. The starburst 

The starburst of NGC 4194 is extended on a scale of 
« 2 kpc (e.g. Wynn- Williams & Becklin 1993; Armus 
et al. 1990) where also 67% of the total ^^CO emission 
can be found. We suggest that overall solid body rota- 
tion of the inner gas (see Fig. 4) provided a dynamical en- 
vironment favourable to a burst of star formation (Keel 
1993). At least 90% of the FIR luminosity can be at- 
tributed to a burst of star formation (Prestwich et al. 
1994; Rowan-Robinson & Crawford 1989) and the high 




2300 2400 '2500 2600 2700 
km/s 



Fig. 12. position velocity cut along the central dust lane, 
from e to a. 



/(60/im)//(100/im) flux ratio (Table 1) indicates that the 
starburst is heating the dust clouds to high temperatures. 

Suchkov et al. 1993 argue that the correlation between 
CO and radio continuum distributions found in galaxies 
is an effect of cosmic ray heating of the molecular gas. 
They suggest that changing values in the ratio of CO and 
radio continuum intensity are indicators of the starburst 
age, with young starbursts having a higher value than 
old bursts since the cosmic rays are more confined to the 
molecular layer if the burst is young. In this model, the 
(low energy) cosmic rays originate in supernove remnants 
(SNRs). Their electron component is responsible for the 
radio continuum emission at cm wavelengths, while their 
proton component heats the gas. This provides a natural 
explanation for the rough correlation between continuum 
and molecular line emission - both are tied to star for- 
mation activity. This relation is, however, also expected 
if the (largely non-thermal) continuum emission simply 
traces SNRs and UV radiation from massive stars plays 
the major role in the heating of the dust and gas. 

It is of course also possible to find molecular gas with- 
out a strong radio continuum counterpart if the gas is 
quiescent and star formation has not (yet) set in. The 
molecular peak e in Fig. 10 may be an example for such a 
region. 

Bonatto et al. (1999) compare NGC 4194 with the 
starburst merger NGC 3256. Their population synthesis 
results suggest that both galaxies have very similar stel- 
lar populations. Starbursts occurcd both at intermediate 
(1-2 Gyr) and young ages. NGC 3256 also has a kpc scale 
starburst (Graham et al. 1984) and is an HR merger with 
relatively faint ^^CO 1-0 emission (e.g. Aalto et al. 1991a) 
and hot dust. 
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4.3. The central gas surface densities and the molecular 

line ratios 

The ^^CO surface brightness of NGC 4194 is lower by an 
order of magnitude compared to other HR mergers such 
as Arp 220 and Mrk 231. Th(; average; gas surface density 
for the extended starburst region of NGC 4194 is 500 — 
1000 Mq pc~^, i.e. more than an order of magnitude below 
the surface density found in ULIRGs for the same CO H2 
conversion factor (e.g. Bryant & Scoville 1999). Address- 
ing the question of the reason for the high TZ posed in 
the introduction, we therefore conclude that the turbu- 
lent, extreme high pressure scenario for compact nuclear 
starburst regions (Aalto et al. 1995) is not well suited for 
explaining the unusual molecular line ratio of NGC 4194. 
Instead, we adopt a similar explanation as for NGC 3256: 
high gas temperatures depopulating the lower levels cou- 
pled with a disrupted and fragmented ISM (Aalto et al. 
1991). One way of testing this scenario is to measure the 
^^CO/^^CO intensity ratio for higher rotational transi- 
tions, and to obtain high resolution images of ^^CO and 
^•^CO like those for the Arp 299 merger. We can also con- 
clude that NGC 4194 is the exception to the general trend 
of increasing central concentration of the molecular cloud 
distribution with increasing /(60/xm)//(100/xm) flux ra- 
tio (Bryant 1997). A study of the intrinsic properties of 
compact and kpc scale HR galaxies is underway (Aalto et 
al., in preparation). 

The star formation rate (SFR) in NGC 4194 has been 
estimated to be; «40 Mq yr^^ (Storchi-Bergmann et al. 
1994), based on Ha emission. This is not much lower than 
the estimate for the HCN-bright ULIRG Arp 220, «110 
Mq yr-i for Arp 220 (Smith et al. 1998). Assuming a sim- 
ilar conversion factor for both galaxies, Arp 220 has about 
an order of magnitude more molecular gas in its central 
region. This would imply a considerably higher star forma- 
tion efficiency in NGC 4194. However, the SFR for Arp 220 
was determined based on its FIR luminosity (see Scoville 
& Soifer 1991 for the method). Thus, the estimates are 
done at different wavelenghts and should be viewed with 
some caution. Indeed, if the SFR of NGC 4194 is deter- 
mined using the FIR luminosity, the result is much lower 
(6-7 Mq yr~^), resulting in star formation efficiencies 
for NGC 4194 and Arp 220 that are comparable. 

The ratio of Lfir/Lco ratio is indeed similar in both 
galaxies implying that both starbursts are equally efficient 
in producing IR photons per unit molecular mass (again 
assuming that CO traces mass in roughly the same way 
in both objects). 

Even if the star formation efficiencies and the 
Lfir/Lco are similar, there may be significant differ- 
ences in the star formation properties of moderate lumi- 
nosity mergers and ULIRGs, since the ISM in the two 
types of galaxies has a different structure. The starburst 
of NGC 4194 takes place in a low pressure environment 
compared to the starbursts at the heart of ULIRGs. The 



lower gas pressure results in a lower mass fraction of dense 
clumps and may allow destructive UV photons to pene- 
trate further into the molecular clouds, where they can 
heat the dust and dissociate the molecules. Then, the 
young, massive stars in NGC 4194 would have a larger im- 
pact on the surrounding gas, disrupting it more effciently. 
This may also account for the relatively high brightness in 
Ha and the hot dust. NGC 4194 has near solar metallic- 
ity, so the apparently efficient heating of the dust cannot 
be explained by a gas-to-dust ratio that is unusually low 
(e.g. Bonatto et al. 1999; Storchi-Bergmann et al. 1994). 

4-4 ■ The Fate of the Molecular Gas 

The elevated ^^CO/^^CO 1-0 molecular line ratio together 

with the IR properties indicate that the molecular gas is 
being strongly affected by the newborn stellar population. 
The estimated star formation rate of NGC 4194 is about 
40 Mq yr""'^ which, if it continues, will have consumed the 
available gas in the center within 40 million years. The 
unusual molecular cloud properties may well be an indi- 
cation that the standard conversion factor is inapplicable, 
so the real gas consumption time scale may be lower. The 
cloud properties will likely be quite similar to those in the 
NGC 3256 model (Aalto et al. 1991a) where the kinetic 
gas temperatures are found to be close to 100 K and the 
clouds are quite small, a few pc in size. Such a molecular 
ISM will radiate more brightly per unit mass than the cool 
disk clouds for which the conversion factor is calibrated. 
Experience from the NGC 3256 model leads us to sug- 
gest that the real molecular mass in the starburst region 
may be a factor of ~5 lower than the mass based on the 
'standard' conversion factor. 

If all of the gas is to be consumed, also the gas at 
higher galactic radii needs to be transported to the center. 
The morphology of the central dust lane with respect to 
the ^^CO peak is very suggestive of gas transport to the 
center. A position-velocity cut (Fig. 12) along the central 
dust lane, from e to a, shows that the gas on the dust lane 
may be moving continuously into the ^^CO peak a. On 
both sides of a along the dust lane the velocities are lower 
by 50 kms~^ than in a, which could imply that gas is 
flowing into a. It is of course also possible that a does not 
lie on the dust lane, but further to the core of the galaxy, 
in which case this may simply be a projection effect. Then 
again, the fact that a is offset from the kinematic center 
argues against this last possibility. There are interesting 
morphological features on the dust lane just east and west 
of a where it seems to flare out into a triangular shape, and 
darken (brighten in the negative Figs.). From Fig. 12 it is 
also clear that at the position of b the emission splits into 
a brighter part at higher velocities and fainter emission at 
lower velocities. 

The final appearance of NGC 4194 may well be that of 
an elliptical with shells and a remaining minor axis dust 
lane. 
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5. Summary 

Wc have obtained a high resolution map of the Medusa 
merger, NGC 4194, in the ^^CO 1-0 line. The main con- 
clusions we draw from this map and single dish ^^CO, 
^^CO and HON observations are as follows: 

1. The molecular distribution in NGC 4194 is extended 
over almost 5kpc. Thus, this moderate luminosity 
merger shows clear morphological differences to typical 
ULIRGs, where an FIR luminosity higher by an order 
of magnitude usually goes along with a very compact 
molecular distribution on a scale of ^ Ikpc. 

2. The structure revealed by our map is complex, consist- 
ing of a rotating central concentration, low level emis- 
sion extending north almost into the base of the tidal 
tail and very prominent arcs and filaments that can be 
identified with two dust lanes visible in the optical. 

3. Along the central dust lane, we can identify several 
clumps that may be Giant Molecular Associations 
(GMAs). Some of these GMAs are positionally co- 
incident with radio continuum features. However the 
strongest molecular concentration visible in "'^^CO is 
offset from the radio continuum peak. This peak may 
represent a second pre-merger nucleus in the system 
or more likely a region of exceptionally strong gas 
concentration and star formation. 

4. NGC 4194 is clearly detected in a single dish ^^CO 
spectrum. Thus, it has a high but not extreme 
12QO/13CO line intensity ratio of ~ 20. The extent of 
the molecular distribution and the relative central sur- 
face density of the gas indicate that the physical mech- 
anism causing the high Tl is different in NGC 4194 
from ULIRGs with high values of TZ. We suggest that a 
scenario with small clouds having high kinetic temper- 
atures may account for the high TZ found in NGC 4194, 
but cannot exclude that diffuse gas may also play a 
role. 

5. The non-detection of the high density tracer HCN in- 
dicates that the amount of dense gas in NGC 4194 
is lower than in ULIRGs. However, the star formation 
efficiency may be nearly the same, with a gas consump- 
tion time of less than 40 million years. Alternatively, 
massive young stars may have a higher impact on the 
fragmented ISM we suggest for NGC 4194. 

6. It is likely that gas is at present flowing along the 
dust lanes toward the nucleus. However, the extended 
molecular distribution is at odds with the classifica- 
tion of NGC 4194 as an advanced major merger. The 
two properties can only be reconciled if the gas flow 
to the center of the system was inefficient compared to 
ULIRGs. Alternatively, NGC 4194 can be interpreted 
as an emerging shell galaxy resulting from a merger 
beween an early type and a spiral galaxy. Comparing 
simulations of such encounters with the morphology of 
NGC 4194 lead us to favour this latter interpretation 
of the Medusa merger history. 



To further our understanding of the merger history of 
NGC 4194, the system needs to be mapped in more molec- 
ular transitions as well as H I and Ha. It represents a class 
of little studied moderate luminosity mergers. More gener- 
ally, knowledge about the parameters that govern the ex- 
tent and efficiency of the starburst resulting from a merger 
and its connection to the progenitors can only be gained 
if the moderate luminosity systems are investigated along 
with the high luminosity ones. If galaxy formation at high 
z happened through mergers, it is essential to understand 
the differences between different types of mergers in the 
local universe, where the systems can be resolved with 
present techniques. 
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